Evidence for the Production of Slow Antiprotonic Hydrogen in Vacuum 
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We present evidence showing how antiprotonic hydrogen, the quasi-stable antiproton (p)-proton 
bound system, has been synthesized following the interaction of antiprotons with the molecular ion 
in a nested Penning trap environment. From a careful analysis of the spatial distributions of 
antiproton annihilation events, evidence is presented for antiprotonic hydrogen production with sub- 
eV kinetic energies in states around n = 70, and with low angular momenta. The slow antiprotonic 
hydrogen may be studied using laser spectroscopic techniques. 
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PACS numbers: 36.10-k, 34.80.Lx, 52.20.Hv 

Studies of the properties of the two-body hydrogenic 
bound states of the stable leptons and baryons have pro- 
duced some of the most precise measurements of physi- 
cal quantities and provided powerful tests of our under- 
standing of the laws of nature. Interest in this area is 
still strong, following the recent production of antihydro- 
gen, H, at low energies [l|, 0]. Accurate comparisons of 
the transitions in hydrogen and antihydrogen are eagerly 
awaited as a stringent test of CPT symmetry. 

Antiprotonic hydrogen (pp) is also of interest. Its level 
structure is similar to that of hydrogen, but with much 
larger binding energies. Precision measurements of its 
spectroscopic properties may allow determination of the 
so-called antiprotonic Rydberg constant and/or the an- 
tiproton/electron mass ratio. 

Although pp has been studied extensively in the past, 
this has exclusively been achieved by stopping antipro- 
tons in liquid or gaseous targets for X-ray spectroscopy 
of inner shell cascades, or for the production of new light 
mesons and baryons (see e.g. [11,13] ). Here we report a 
radically new method of pp production resulting in emis- 
sion almost at rest in vacuum. This has been achieved 
using a chemical reaction between antiprotons and molec- 
ular hydrogen ions (H+) in the ATHENA Penning trap 
apparatus [3, Q . This advance has opened the way for 
laser spectroscopic studies of pp, or other antiprotonic 



systems formed by p interactions with HD + or , akin 
to those successfully deployed in the study of antiprotonic 
helium (see e.g. 0,0])- 

The experiments were made possible by the availabil- 
ity of a high-quality low energy p beam delivered by 
the CERN Antiproton Decelerator to the ATHENA H 
apparatus. The latter contained a multi-electrode sys- 
tem of cylindrical Penning traps, 2.5 cm in diameter and 
~ 90 cm in length immersed in an axial magnetic field 
of 3 T. The residual pressure of - 10" 12 Torr, in the 15 
K cryogenic environment of the trap, was due to hydro- 
gen and helium gases. The central region contained the 
mixing trap: a nested Penning trap, approximately 10 
cm long, that allowed positrons, e + , and antiprotons to 
be confined simultaneously. For H production the mix- 
ing trap contained a spheroidal cloud of ~ 3.5 x 10 7 e + . 
Around 10 4 antiprotons were injected into this plasma, 
with the resulting p annihilations monitored for 60 s by 
position sensitive detectors [5j, [8( . These registered the 
passage of the charged pions, to localize annihilation ver- 
tices which were due, not only to H formation followed 
by annihilation on the electrode surface 0, Q , but also 
p a nnihilation following transport to the electrode walls 
[1 01 ] , and annihilation following interactions with residual 
as atoms or ions present in the trap. It was shown in 
that the vertex data were predominantly H anni- 
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Figure 1: r — z scatter plot and radial densities (^4tO of 
the annihilation vertices for: (a) hot mixing; (b) cold mix- 
ing. The dashed black line indicates the position of the trap 
wall; the red semi-ellipse shows the section of the e + plasma. 
The green radial densities are for the corresponding central 
z-region events (inside the green lines) whilst the data in blue 
are for the 2 lateral z-regions, normalized for r > 1.25 cm. 



hilations during so-called "cold mixing" (CM), when the 
e + cloud was held at the trap ambient of 15 K. In con- 
trast for "hot mixing" (HM), when the e + were heated 
to a temperature, T e , of several thousand K (here 8000 
K) [ll|, |T3| , H formation was suppressed and the p anni- 
hilations were mainly a result of collisions with trapped 
positive ions. It is this effect (which is also present for 
CM) that is addressed here. 

FiglH shows r — z scatter plots for annihilation vertices 
taken under HM and CM conditions. Here the radial 
positions, r, (i.e. the distance from the trap axis) of 
the events are plotted versus their axial coordinates, z. 
Also shown are the attendant radial density distributions 
(r cF )" ^ ne distributions are broadened by the uncer- 
tainties in the vertex determination (around 1.8 mm in 
the z-direction and 3.5 mm in the transverse dimensions) 
caused mainly by the inability of the ATHENA detector 
to reconstruct the curved trajectories of the pions in the 
3 T magnetic field. The present HM results are for the 
highest T e achieved by ATHENA. 

There are striking differences between the two r — z 
plots. Besides the H annihilations on the trap wall cen- 
tred around r = 1.25 cm (CM only), there are events lo- 
calized at smaller radii which dominate in HM (FigJTJi), 
but which are also evident in CM (FigJTJ)). Examin- 
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Figure 2: Experimental radial distribution of p annihila- 
tion vertices (black histogram) for HM ( — 1.5 cm< z < 
1.5 cm) with a Monte Carlo simulation (T e = 8000 K, Vth 
= 5600 ms~ , generation on the surface of a spheroid with 
z p — 16 mm and r p = 1 mm rotating with a frequency of 



300 kHz, i.e. 



Vta 



2000 ms" 



see text for details. Re- 



sults of simulations with different mean lifetimes are shown: 
green, r = 0.8/is (xled = 2.78); red, r = 1.1 /is {xled = 1-48); 



blue, 



1.4*18 {xled = 2.14). 



ing FigQJ), the radial density distributions for CM at 
0.5 cm) behave quite differently be- 



small radii 



< 



tween the central {\z\ < 0.5 cm) and its adjacent regions 
(0.5 cm < |z| < 1.5 cm). Moreover, the shape of the dis- 
tribution for the events in the central region resembles 
that for HM. 

In FigHand FigsKi,b the radial (^) and axial (^) 
annihilation distributions are plotted for the HM case. 
FiglSJ: shows the axial distribution for CM for events with 
r < 0.5 cm being notably narrower than the HM case in 
FigEb. 

In order to determine the characteristics of the near- 
axis events the possibility that they are due to H has 
been investigated. In ATHENA, H was detected by the 
coincidence in space and in time (within 5 /is) of p and 
e + annihilations. This was achieved for events having a 
charged pion vertex accompaincd by a pair of 511 keV 
photons with the angle between them denoted # 77 [j], 0] . 
Examination of the distributions of cos # 77 indicates that 
the trap centre events are not due to H, except for a small 
fraction in CM where the long tails in Fig|3fc are due to 
poor reconstructions of H annihilations on the trap wall. 
No H annihilations occur for HM. Moreover, cos# 77 dis- 
tributions for CM in the three z-regions in FigQJ) suggest 
that the fraction of non-H annihilations on the wall in the 
central z-region is insignificant. 

The distributions in Figs 03 [2] also show that in HM the 
annihilation distribution is extended and has a conspic- 
uous fraction out to the trap wall. We interpret these 
features as evidence for the production of slow antipro- 
tonic hydrogen as follows: 

i) the limited axial range of the vertices, when com- 
pared to the total length of the nested trap, indi- 
cates that p do not annihilate in-flight on residual 
gas which is present throughout; 



3 





80 


=3 


60 


o 
O 


40 




20 







<n 


80 


c 

=3 


60 


o 
O 


40 




20 







tn 


200 


c 

=3 


150 


o 
O 


100 




50 








r > 1 .0 cm ^-rfrffH- 




n (a) 


£vr-rT7T! . i .... i .... i .... i .. ..■ . 


, , , i , , , , i , .TTTrW-?: 



; r<0.5cm 

L=^j=4i=F | T 1 -rT i .... i .... i .. . 


i tfWfl (b) 

. i .... i .... i .... 1 7r?=ui n, n. 




'_ r<0.5cm 


f\ (c) 



2 

z (cm) 



Figure 3: Experimental axial distributions: hot mixing for 
events near the trap wall (a), and near the trap axis (b) and 
for cold mixing for events near the trap axis (c). In (a,b): the 
red line is the simulation with the parameters of Fig[2] and 
with lifetime r = 1.1/is. In (c) the red line is the simulation 
with the parameters of Fig(4^i. 



ii) the events cannot correspond to in-flight annihila- 
tion on trapped positive ions since the latter are 
only present near the e + cloud whilst the annihila- 
tion events are radially diffuse; 

iii) positive ions can capture p in the central part of 
the recombination trap. In the case of pHe + the 
residual electron would be rapidly ejected, in the 
majority of cases in less than 10 ns [13j, leaving a 
charged system that would annihilate very near its 
point of formation. Since this is inconsistent with 
our observations, p capture by helium ions can be 
excluded. 

Further information on the inferred antiprotonic sys- 
tem formed from capture by a positive ion can be ob- 
tained by exploiting the different charged pion multiplic- 
ities expected for p annihilation on a proton or neutron. 
Tab|T] shows the ratios, R23, of the number of the recon- 
structed annihilation vertices having two charged pion 
tracks to those with three tracks for different data sam- 
ples. 



Data set 


Ratio f?23 on wall 


Ratio -R23 at centre 


Cold mixing 


1.35±0.01 


1.22±0.04 


Hot mixing 


1.38±0.10 


1.17±0.04 


antiprotons only 


1.40±0.03 




Monte Carlo pp 


1.19 ±0.01 


1.19 ±0.01 



Table I: Experimental and Monte Carlo results for the number 
of charged pion tracks due to p annihilations. 



From TabU the R23 values for annihilation on the trap 
wall for all the samples agree, within uncertainties, but 
differ from those for the trap centre by 4 standard de- 
viations. This is not due to geometry, since the Monte 
Carlo simulations assuming the pp system give the same 
result both for the "wall" and "centre" annihilations. It 
is likely that the trap centre events are due to pp annihi- 
lation, since the pp Monte Carlo result agrees well with 
the experiment. This, together with the three constraints 
described above, suggest that pp is responsible for the ob- 
served annihilation distributions. Thus, the data can be 
further analysed to search for consistency with the con- 
ditions pertaining to the production region at the two 
different positron temperatures, the pp lifetime and fi- 
nally the energetics of the formation reaction. 

First, pp is not confined by the electromagnetic fields 
of the traps and if formed in a metastable state 0, EH, 
, 1^1 it can decay in flight far from its point of forma- 
tion, perhaps even on the trap wall. The convolution of 
the pp lifetimes and velocities governs the observed anni- 
hilation distributions. Moreover, Figs. [3Jd and c clearly 
indicate that the temperature of the e + cloud influences 
the spatial origin of the pp. 

We have attempted to generate the initial position of 
the pp using a Monte Carlo simulation which takes into 
account the radius, r p , and the axial half-length, z p , of 
the e + spheroid, both deduced by means of the non- 
destructive technique described in [ill . [T3 | . Parameters 
obtained using this technique allow the plasma rotation 
frequency to be extracted. 

In this simulation, the spheroid was characterized by 
r p = 1 mm and z p = 16 mm. For the CM case the pp 
was generated in a region with a fixed radial position at 
r = r p = 1 mm and with a Gaussian distribution along 
the axis centred at the symmetry plane of the plasma 
with cr = 2.5 mm. This gave the best fit to the data. For 
the HM case, however, pp was generated with a = 10 mm, 
though limited to the length of e + plasma. It is notable 
that, for the HM case, the simulated annihilation distri- 
butions were not strongly dependent upon the assumed 
starting conditions, taking into account our experimental 
resolution. 

The velocity of the pp was generated from the sum of 
a thermal Maxwellian distribution, Vth, smd the tangen- 
tial velocity, Vtang, induced by the ExB plasma rotation 
as vtang = E x B/|B| 2 . Following this prescription, the 
mean radial kinetic energy of the pp is about 40 meV for 
CM (15 K), and dominated by the effect of the plasma ro- 
tation, and about 700 meV in the HM case (8000 K), and 
dominated by the plasma temperature. An exponential 
decay law for the pp lifetime distribution was assumed 
such that its mean lifetime was determined by fitting the 
simulations to the observed data. 

The simulated radial and axial annihilation distribu- 
tions are plotted with the HM data in Fig^and Fig[3^,b. 
The agreement is good and the best fit was obtained with 
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Figure 4: Experimental radial distribution for cold mixing 
2003 (a) and cold mixing 2002 (b), obtained by subtracting 
the H contribution (see text). The red lines are the simula- 
tion results. In (a): T=15 K, vth = 250 ms' 1 , generation 
from a spheroid with z v = 16 mm and r p = 1 mm rotating 
with a frequency of 300 kHz, i.e. Vtang = 2000 ms" 1 ; mean 
lifetime is r= 1.1 /is. In (b): same parameters as in (a) ex- 
cept r p = 2.5 mm and the rotation frequency is 80 kHz, i.e. 
vtang = 1300 ms _1 . In (b) the green line corresponding to 
the parameters of the red line in (a) is shown for comparison. 



a mean lifetime of (1.1 ± 0.1)/j,s. The sensitivity of the fit 
to this lifetime is illustrated in Fig[3]by the clear discord 
between the experimental data and the simulations when 
lifetimes of 0.8 /is and 1.4 /its were used in the latter. The 
fits imply that about 25% of the antiprotonic hydrogen 
atoms annihilate on the trap wall in HM. 

To isolate the radial distribution of the non-fi anni- 
hilations in the central ^-region for the CM sample the 
normalized distribution, as evaluated from Fig[TjD for the 
two lateral z-regions, was subtracted from the central 
one. Since, as noted above, non-H annihilations on the 
wall in the central z-region are insignificant for CM, the 
radial distributions were normalised for r > 1.25 cm. 

The results are plotted in Figf4^. The Monte Carlo 
simulated events, assuming the same lifetime extracted 
from fits to the HM data, also show good agreement for 
CM. In this case, the simulation indicates that < 0.5% 
of the pp reaches the trap surface. 

A further test of our model has been obtained by exam- 
ining a sample of CM data acquired by ATHENA in 2002 
with a e + plasma radius of 2.5 mm (and hence a differ- 
ent rotation frequency). Figj4j3 shows the radial vertex 
distribution of this sample together with Monte Carlo 
simulations for two values of the radius of the pp source. 
The simulated events for the 2.5 mm source are in much 
better agreement with the experimental data than those 
for 1.0 mm. This is illustrative of the sensitivity of the 
model to the spatial origin of the pp for CM samples. 

A possible explanation of the experimental differences 
between the pp distributions in CM and HM lies in the 
nature of the thermal equilibrium state of the combina- 
tion of a e + plasma with an admixture of ions. The 
physics of the radial separation of the different species in 
two-component plasmas has been given elsewhere 171.118] 



and has been experimentally observed for a mixture of 
positrons and 9 Be + in [l9j . For our experimental condi- 
tions, assuming thermal equilibrium, the centrifugal po- 
tential barrier is of the order of 10 meV. Thus, for CM at 
15 K, the thermal energy of the ions means that they will 
be partially separated from the e + and confined near the 
equatorial region of the plasma. However, at a e + tem- 
perature of 8000 K the barrier is negligible and the ions 
will be present throughout the plasma. 

We note here that the experimental data cannot be re- 
produced by simulation if it is assumed that the pp gains 
a recoil energy of the order of 1 eV or higher. Thus, it 
is contended that pp is being produced in a recoil-free 
collision of an antiproton with a positive ion such that 
its velocity is determined predominantly by the thermal 
and plasma environment. This insight, together with the 
constraints i)-iii) noted above, suggest that the only pos- 
sible collision partner for the p is the molecular ion, H^~. 
Thus, the inferred pp production mechanism is, 



P 



H 2 + 



pp{n, I) + H. 



(1) 



Hence we believe that ATHENA has observed around 
100 antiprotonic hydrogen annihilations every 60 s p in- 
jection cycle for the CM and HM conditions. In an exper- 
iment in which the number of ions present in our nested 
well was counted by measuring the charge collected fol- 
lowing emptying of the trap it has been deduced that the 
trap contained Hj ions. These probably arose as a re- 
sult of the positron loading procedure |2(j, in which the 
positrons could collide with H 2 residual gas as they were 
slowly squeezed into the mixing region. Ions may also 
be produced and trapped during p loading (see also [211 ]) 



and we estimate that around 10 4 — 10 5 ions were present 
under typical ambient conditions. It is straightforward to 
show that this ion density, together with the observed pp 
production rate and the p speeds used in the simulation 
are consistent with calculated cross sections for reaction 
(1) 0. 

Agreement with [16[ does not, unfortunately, extend 
to the most likely principal quantum numbers of the an- 
tiprotonic h ydr ogen atoms produced in the reaction. The 
calculation [161 ] finds production peaked around n = 34, 
in the presence of substantial pp recoil. The latter is 
contrary to observation. Simple kinematics relating to 
near zero-energy p — H^ collisions suggest that n = 68 
should dominate, with the liberated hydrogen atom in its 
ground state. In this case the lifetime of 1.1 /is extracted 
from the simulations implies that production in low an- 
gular momentum states (I < 10) is favoured, since the 
radiative lifetime to an I — state (which is followed by 
prompt annihilation) weighted by the statistical (2/ + 1 ) 
distribution is around 15 fj,s for n = 68. The dominance 
of low I is intuitive for such a slow collision in which the 
molecular ion will be severely polarised by the incoming 
p, resulting in an almost collinear collision system. 
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In conclusion we have presented evidence for the pro- 
duction of antiprotonic hydrogen, in vacuum, with sub- 
eV kinetic energies and in a metastable state. Given the 
capability of accumulating 10 s ions in tens of sec- 
onds and storing ~5x 10 6 antiprotons in some minutes 
22l 23[, our result opens up the possibility of perform- 
ing detailed spectroscopic measurements on antiprotonic 
hydrogen as a probe of fundamental constants and sym- 
metries. 
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